Optical sensor array system and method for parallel processing of chemical and biochemical information

ABSTRACT

Total analysis systems and methods for simultaneously monitoring a suite of biological and/or chemical species in water and/or other process systems are disclosed. The system provides a sample-volume controlled sensor array comprising a fluid delivery device and a plurality of optical sensor elements for determining the presence and total concentrations of multiple analytes in the process system simultaneously. Image identification algorithms are provided for identifying the analytes based on image intensity, color pattern, positional arrangement, and the like. The methods incorporate multivariate optimization algorithms to analyze multiple sensor responses. This produces analytical results that are typically difficult to obtain without full system or variable compensation. The improved array response may then be utilized to measure, monitor, and control the concentration of analytes in the chemical or biological sample or water system.

FIELD OF INVENTION

The present invention relates generally to a chemical sensor array, and more particularly relates to a new and improved system, method and apparatus for delivering a fluid sample to a chemical sensor array, and parallel processing chemical and biochemical information from a multiplicity of sensor elements of the sensor array.

BACKGROUND OF THE INVENTION

Many chemical and biological measurements need to be performed in locations outside of a fully equipped analytical facility. This requires systems that are portable and miniaturized so that they can be transported to locations where rapid test response is required for process or water quality monitoring, or they can be deployed in a medical environment to provide rapid test results for certain biological or biochemical species of interest. These chemical and biological analyses can be performed individually using single tests with post-test optimization to enhance the quality or accuracy of the results, but this serial approach has inherent flaws because multidimensional interactions are difficult to fully compensate using a serial approach. Additionally, this approach can be time consuming and can produce erroneous results. Introducing operator or system errors when tests are performed on different platforms or at different times further complicates this system. The best way to overcome this limitation is to perform all the desired measurements simultaneously on the same platform, but current state of the art does not provide a fully integrated platform for such measurements.

Current electrochemical array technology allows operators to perform an arrayed test at a single time, but this is limited to those materials that respond to an electrochemical stimulus. This usually involves measurement techniques like anodic stripping voltametry or cyclic voltametry methods, or incorporates chemically responsive material into an electrochemical detector, e.g., an ion-specific electrode (ISE). This system, although productive for some systems is limited by many of the common limitations of electrochemical systems, e.g., systematic problems at low and high ionic strength that effect electrochemical potentials. Additionally, some of these systems can suffer from serious cross-reactivities or interferences, e.g., the cross-reactivity of common oxyanions or small cations like sodium, lithium, and potassium.

There are other test platforms that could provide small-scale array measurements that are based on optical or spectral measurements. These can be optical detection from multi-flow wet chemical analysis, or they can be portable versions of classic laboratory measures, e.g., portable atomic absorption spectrometer units. These systems are often limited by mechanics required for fluid flow and maintenance, or from cumbersome equipment like portable atomic absorption spectrometer systems that, although theoretically transportable, have proven less mobile in practice. There is also mention of miniaturizing additional lab systems like Inductively Coupled Plasma—Atomic Emission Spectroscopy or Mass Spectroscopy, but these methods are difficult to adapt into portable handheld, or field deployable systems.

A proposed alternative is to use an optical platform based on well-characterized, chemical responses of optical sensor films. Such a system uses solid, chemically responsive films that respond to analyte concentrations by changing their absorbance values at an optimized wavelength. This platform can be extended to incorporate sensor test elements for all known interfering or cross-reactive species for a particular test matrix, as well as account for test limitations at the extremes of the test sample conditions, e.g., high and low ionic strength as well as high and low buffer strengths. This system has the added benefit of providing a small test platform that can be formed into an array specifically designed to measure test elements that require specific deconvolution analysis.

Optical chemical sensors also fall into two general classes, reversible and irreversible. Fully reversible sensors equilibrate rapidly to the concentration of the target analyte in the test fluid and their response changes when the analyte concentration changes. Examples of reversible sensors are polymer film pH sensors and ion selective optodes (ISO). In contrast, an irreversible sensor will continue to respond to the analyte in the test fluid until the responsive reagent in the sensor has been exhausted, i.e., the total amount of analyte available to the sensor rather than the analyte concentration in the sample. Many non-ISO type sensors belong to this category.

Since the reagent in the reversible sensors is at chemical equilibrium with the analyte in the sample, the exposure of the sensor film to the sample alters the analyte concentration if the sample volume is finite. This requires that the reversible sensor films be exposed to either a large excess amount of sample volume or a given amount of sample volume. In the latter case, a correction can be made to reduce errors due to the finite volume effect. Similarly, the irreversible sensor requires sample volume control so that the sensor response reflects the analyte concentration in the controlled volume of the test fluid.

A sensor array designed for quantitative analysis may not yield satisfactory results by just immersing the array element into a liquid sample because the above mention reasons. For a sensor array consisting of both reversible and irreversible sensors, the sample volume that each sensor region is exposed to has to be controlled. Moreover, volume regulation also helps prevent sensor-to-sensor cross contamination. In this invention, sensor film compositions are designed to be at their optimized performance when they are exposed to fixed sample volumes.

Optical sensor arrays that are composed of irreversible sensors, or a combination of irreversible and reversible sensors, must have some form of fluidic control that delivers a controlled volume of test fluid to each sensor element. Most systems available today use some form of pump or mechanical multi-addition system to deliver these controlled volumes, e.g., robotic addition to multi-well plates. These systems are often cumbersome and require mechanical and electrical components that are rarely field robust and suitable for remote testing in harsh environments.

In recent years, capillary effect has been exploited for fluidic designs. One drawback associated with this passive mechanism is that it relies on the use of absorbent or wicking materials. This makes it difficult to fabricate a device to deliver a small amount of sample to a large number of locations. Instead, the absorbent material is an integrated part of sensing or reaction matrix. Liquid delivered to the site result in only wetting the materials within the matrix. For sensor array applications, dosing a given amount of liquid sample to multiple locations is more desirable.

Even though a large number of publications and patents have been devoted to the development of sensor methods, reagents, and equipment to replace the traditional wet chemistry methods, a need remains for an economical and convenient field deployable sensor system for simultaneous detection of multiple analytes.

SUMMARY OF THE INVENTION

A sample-volume controlled sensor system for simultaneously measuring multiple analyte concentrations in chemical or biological substances such as water systems comprising a set of reversible and irreversible analyte-responsive sensor elements which are selected to change at least one optical property in response to chemical, biological, or environmental stimuli, including at least one reference region that serves as an internal optical and position standard, and a light source for directing light onto an array of sensor elements. A detector-based imaging device is provided which can adjust its performance to the position and spectral profile of the light source, and then converts this imaged response into a digital record. Image identification algorithms are provided to identify the test composition on the element by one of many configurations based on image intensity, color pattern, arrangement, and the like. A software-based optimization algorithm is provided which incorporates responses from the sensor array and produces optimized results unavailable without full system and variable compensation.

In this invention, a device comprising channels and reservoirs capable of delivering a controlled amount of a liquid sample to multiple reservoirs containing an array of sensor elements within a given time period is also disclosed. The driving force for liquid transport within the device is mainly capillary force generated by the surface energy of the liquid and channel/reservoir wall interface. This device does not rely on the use of any wicking materials and can be produced inexpensively using readily available materials.

Also disclosed is a total analysis method for monitoring a suite of biological and chemical species in water and process systems. The system provides a sample-volume controlled array of reversible and irreversible optical sensors to determine total concentrations of multiple analytes simultaneously. The method involves exposing the sensor array to a multiple analyte-containing media, and recording the sensor array response as a digital record. The sensor response is processed to reduce noise and interference, and multivariate analysis is applied to improve the array response. The improved array response is then utilized to measure, monitor, and control the concentration of analytes in the chemical or biological substance or water system.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a perspective view of an assembly stack for the capillary-flow fluid delivery sampling device in accordance with an exemplary embodiment of the present invention;

FIG. 2 illustrates a branched series channel-reservoir configuration in accordance with another exemplary embodiment of the present invention;

FIG. 3 illustrates a parallel-series channel-reservoir configuration in accordance with another exemplary embodiment of the present invention;

FIG. 4 illustrates a fluid delivery configuration with a waste reservoir and a delay channel in accordance with another exemplary embodiment of the present invention;

FIG. 5 is an image of chlorine sensitive films from right to left, chlorine concentrations are 1, 2, 4, 5, 10, and 50 ppm;

FIG. 6 illustrates the calibration curve for chlorine determination;

FIG. 7 is an image of alkalinity sensitive films;

FIG. 8 is a graph that illustrates [(R_(w)−R)² +(G_(w)−G)²+B²]^(1/2)/B_(w)calculated from a digital image captured with a scanner plotted a function of solution alkalinity;

FIG. 9 is a graph that illustrates the correlation of [(Rw−R)²+(Gw−G)²+B²]^(1/2)/B_(w) to absorbance measured with a optical reflectance probe at 650 nm;

FIG. 10 is a graph that illustrates [(R_(w)−R)²+(G_(w)−G)²+B²]^(1/2)/B_(w)calculated from a digital image captured by a color digital camera plotted a function of solution alkalinity;

FIG. 11 is a graph that illustrates performance of the camera vs. scanner;

FIG. 12 is a graph that illustrates the multivariate calibration of pH determination;

FIG. 13 shows a pH sensor film with defects;

FIG. 14 is a graph that illustrates the effects of film defect on RGB values;

FIG. 15 is a graph that illustrates the rejection of a group of pixel by standard deviation criteria;

FIG. 16 is a graph that illustrates the calibration curve for Ca sensor film;

FIG. 17 is a graph that illustrates kinetic response of molybdate sensor film at different temperatures;

FIG. 18 comprises digital images showing color change of the molybdate sensor film during the course of its exposure to a water sample at 25.3° C.;

FIG. 19 is a graph that illustrates temperature effect on molybdate sensor response;

FIG. 20 shows the channel and reservoir layout for Example 1;

FIG. 21 shows comparison of predicted fill time with experimental fill time data for Example 1;

FIG. 22 is a graph that shows average fill time for 54 μl reservoirs with error bars obtained from six tests in Example 2; and

FIG. 23 is a calibration curve for the magnesium sensitive film obtained with a sample-volume controlled sampling device detailed in Example 9.

DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS

This invention describes a test array platform that incorporates noise reduction, interference reduction, improvements from multi-chemistry response, multivariate analysis and a flexible platform, that allows test arrays to be custom designed to provide optimized response with minimized systematic or operator errors. This array-based testing platform is based on chemically responsive optical sensor elements that can be incorporated into rugged, field deployable systems composed of a simple detection array, and these systems can be easily interfaced with computer or electronic units that can perform the complex analysis required to provide optimized measurements at non-laboratory locations.

One aspect of the present invention is the recognition that optical sensor films can be developed to account for most systematic variables that effect chemical and biological analyte measurements. These chemical and biological sensors systems contain several functional components. One component is a sensor material that responds to a change in the environment. Examples of such sensor materials are analyte-responsive polymers, biomembranes, sol-gels, and some others. In the case of optical sensors, the sensor material should maintain adequate optical transparency or loss of transparency for the chemically responsive composite to be monitored using optical transmission, reflection, dispersion, fluorescence, or any other common optical method known in the art. Another component is the electronic system that provides means of measuring the change in the sensor material upon environmental exposure. Interactions of the environment with the material are converted into an analytically useful signal using an appropriate transduction mechanism such as optical sensing. This arrayed optical detection platform is interfaced with a “smart system” that compensates for interferences, environmental variations, etc., and performs noise reduction and test optimizations and produces final results of greater quality than can be obtained with a system that is not as fully integrated.

The following sections describe the components of this total analysis system, with examples of how each component can provide increased improvement when applied to a complete array test platform. The total analysis system is a product of combining each of these improved elements to produce a single system with enhanced performance that results from the simultaneous combination of the individual elements.

Sensor Array

An optical sensor array has a set of analyte-responsive elements where the sensor elements respond to analyte concentrations by changing color or other optical properties upon exposure to a sample. The number of total sensor elements and the type of sensor elements can be selected to meet a specific system analysis need. As a non-limiting example, one sensor array for water analysis comprises optical sensor elements responding to the following analytes: Alkalinity, pH, chlorine, hardness, sulfite, and phosphate.

Suitable sensor types for use in the present invention are described in our co-pending patent applications entitled “Material Compositions for Sensors for Determination of Chemical Species at Trace Concentrations and Method of Using Sensors” and “Self-Contained Phosphate Sensors and Method for Using Same” filed on the same date as the present application, the disclosures of which are hereby incorporated by reference herein.

Optical sensor arrays comprise an array of solid sensor elements deposited on a solid substrate. The solid element can contain a single or multiple components. One or all components in the solid sensor element could be water-soluble. Combinations of different solubility characteristics of the components in the sensor element can be chosen to enhance sensor array performance. As an example, the sensor element can be prepared from a hydrogel polymer containing water soluble reagents that respond to analyte concentration.

A binder may be used to enhance adhesion of the solid element to the substrate. A liquid spreading material, e.g., a surfactant, may be added to the solid element to improve the wetting properties of the sensor regions. The liquid spreading material may be placed in between the solid element and the substrate, or in other configurations such as on the top of or surrounding the element. Commonly practiced methods that are suitable to prepare the sensor array needed for this invention include the methods for manufacturing indicator paper strips and polymer film sensors as described by Zolotov et al. in “Chemical Test Methods of Analysis” in Wilson & Wilson's Comprehensive Analytical Chemistry, 2002, the disclosure of which is hereby incorporated by reference herein.

Array Light Source—Detector Combinations

There are many light source/detector combinations suitable to measure sensor responses on an optical array. For example, our prior U.S. patent application Ser. No. 10/760,438 filed Jan. 20, 2004 describes a handheld device with a disposable element for chemical analysis of multiple analytes, the disclosure of which is also hereby incorporated by reference herein.

Turning to our present invention, the invention pertains to novel systems and methods for simultaneously detecting parallel film responses from a plurality of sensor elements. The table below shows sources for UV-visible-near-IR ranges for the applications in conjunction with optical sensor array system, and for parallel processing of chemical and biochemical information suitable for use in the present invention. It is understood that other less conventional light sources that emit radiation in the spectral range of interest such as sun, organic light emitting diodes, indoor room lights, products of bioluminescence reaction, emission of electronic equipment such as computer monitors, PDA monitors, displays of cell phones, pagers, radioluminescent sources, and any other light sources known or later developed in the art could also be used without departing from the scope of the present invention.

Light Sources Useful For Optical Sensors Source Spectral range of emission (nm) Continuous wave sources: Xenon arc lamp 200-1000 Mercury arc lamp 250-600 Deuterium lamp 180-420 Tungsten lamp 320-2500 Light emitting diodes different diodes cover range from about 250 to 1500 nm Diode lasers different diode lasers cover range from about 400 to 1500 nm Argon ion laser several lines over 350-514 nm Helium-neon laser several lines over 543-633 nm Krypton laser several lines over 530-676 nm Pulsed sources: Nitrogen laser 337 nm Nd: YAG laser fundamental - 1064, frequency doubled - 532, tripled - 355 Ti: Sapphire laser 720-1000, frequency doubled 360-500 Dye lasers 360-990 frequency doubled 235 to 345

Possible detectors include vacuum or solid state and single or multichannel detectors. Vacuum detectors are phototubes and photomultiplier tubes (PMT). Solid-state detectors include photodiodes, photodiode arrays, charge-coupled devices (CCDs), charge-injection devices (CIDs), and avalanche photodiodes. Multichannel detectors include arrays of individual detectors such as photodiode arrays, PMT arrays. Also, CCDs, CIDs, CMOS, and other types of multichannel detectors are available. Each element has its intrinsic advantages and disadvantages and can be combined to produce a light source detector platform suited to the particular need in a specific application. Similarly, it is possible to combine more than one light source or detector to monitor different types of responses on sensor films, and then combine them into a common array platform in a manner known in the art.

As an example, a color image of a subject can be recorded by illuminating with any of the cited or envisioned light sources and captured by a digital scanner or camera. A CCD color sensor in the digital camera measures the intensities of the three primary colors (red, green, and blue) of the subject. The Red-Green-Blue (RGB) color intensity value of each pixel is recorded in a digital file. The color depth, or the range of RGB values, is usually 0 to 255. Colors can also be measured by the CCD color sensor in a digital scanner using a white light source for illumination and some form of simple color light detector. Some digital scanners, however, use three LEDs (red, green, and blue) to irradiate the subject for the color measurement. Unlike digital cameras, most digital scanners provide a 48-bit, or greater, color resolution. In this color model, the color of each pixel is quantified by RGB values in the range of 0 to 65025. The spectral ranges of the three primary colors measured by the digital camera and scanner vary slightly from model to model. The spectral response of a typical CCD color sensor is 460±40 nm, 540±40 nm, and 645+/−55 nm, respectively. We use a digital imaging device for quantitative colorimetric analysis. In this aspect, the digital imaging device is equivalent to multiple, three-color LED/photodiode pairs.

Fluidic Delivery System

FIG. 1 illustrates a fluidic delivery device 10 in accordance with an exemplary embodiment of the present invention. The delivery device 10 transports a controlled amount of a liquid sample, in metered quantities, to multiple reservoirs 8 in order to effect a chemical reaction between the sample fluid and the sensor elements (not shown) connected to the reservoirs 8. As shown in FIG. 1, the fluid delivery device 10 comprises a top cover layer 2, a middle channel layer 4, a bottom sampler-substrate binding (i.e. gasket) layer 6, a fluid entry port 12, and an associated plastic entry port wall ring 11. A plurality of grooves or channels 5 are formed on the channel layer 4 for directing the sample fluid from the fluid entry port 12 to the reservoirs 8. A plurality of channels is formed when the cover layer 2 is bound to the channel layer 4. A series of vent holes 7 are added to assure complete fluid flow through the channel system.

Many commercially available hydrophilic films can be chosen as the top layer to fabricate the device disclosed in this invention. Some films have a heat sealable adhesive deposited on the hydrophilic side. For those films without adhesives, standard-bonding methods can be used to laminate the cover layer to the channel layer such as ultrasonic welding and adhesive transfer bonding. Hydrophilic films may be both heat sealable and pressure sensitive adhesive.

The channel layer 4 can be fabricated using conventional plastic processing methods such as injection molding, hot embossing, and micro machining. Many plastic materials that have water contact angle in the range of 40 to 85 degree can be used for the channel layer. For example, polycarbonate and acrylic are suitable materials for this application.

The sampler binding layer 6 can be any material that has durometer number around 40 Shore A and provides seal to a flat substrate by surface wetting, conformal contact and/or adhesive bonding. Non-limiting exemplary materials for this application are silicone and synthetic rubbers, and thermoplastic elastomers. The sampler binding layer 6 could be an extruded silicone sheet. A double side adhesive material could also be used. The sampler binding layer 6 can be bonded to the channel layer using adhesives, although by choosing a high-heat plastic material, such as polycarbonate or Ultem, insert or two-part molding of the sampler binding material to the channel layer can be done.

The cover layer 2 provides a super hydrophilic surface for the fluidic channel, which contributes largely to the overall capillary force that drives the fluid to flow through the channel. The cover layer also accommodates a plurality of small holes 7, one over each reservoir 8, to allow for the passage of air, which is displaced out from the reservoir by the incoming liquid. Due to the capillary force driving the fluid through the channels 5, no pumps and valves are required to deliver a given amount of liquid sample or reagent from the sample entry port 12 to the multiple reservoirs 8 within a predefined sequence. As a result, the device 10 can be efficiently fabricated by inexpensive plastic processing methods. The fluid delivery device 10 may then be integrated as a component of a chemical or biological sensor array system to transport and dose, in metered quantities, a sample liquid into the reservoirs 8 in order to complete a reaction with the associated sensor elements.

In order to effectively fill a reservoir 8 from a channel 5 that has a smaller capillary dimension than the reservoir itself is not a trivial matter. The transition zone from the channels to the reservoir is likely to act like a capillary barrier in obstructing the passage of liquid from the end of the channel to the reservoir. An external force, such as gravitational force, may be needed to overcome this barrier. In other cases, channel and reservoir parameters are carefully balanced to shorten the transition time and yet avoid overflow through the vent holes 7.

The present invention overcomes the capillary barrier transition problem described above by implementing the following design features. First, a one-sided super hydrophilic film is chosen for the cover layer, which enables the liquid to advantageously accumulate over the entire top wall of the reservoir to form a pendant droplet. As the droplet grows, gravitational force enables it to reach the bottom wall of the reservoir and then capillary force originated from the four walls of the reservoir drives the liquid to fill the entire reservoir. As the reservoir completely fills, the vent hole serves as a capillary barrier, impeding the liquid from flowing through it to reach the top surface of the cover layer, which is designed to be hydrophobic. Secondly, we adjust the fluid flow resistance through optimization of channel and reservoir geometric parameters, hydrostatic pressure from the entry port, and capillary pressure of the channel in order to achieve a desirable fill volume and filling sequence. In addition, the sampler binding layer 6, which appears as part of the sidewalls of the reservoir 8, creates a capillary barrier as the pendant drop grows. Therefore, a careful design of the channel and reservoir parameters, including the sampler binding thickness, is important to overcome this barrier to ensure complete filling of each reservoir.

Referring again to FIG. 1, each channel 5 can feed single or multiple reservoirs 8. If a channel is required to feed multiple reservoirs due to space consideration, a simple branched-type structure as shown in FIG. 1 may be constructed to help prevent trapping of air bubbles in the channel 5. In the branched-type configuration, the filling sequence of reservoirs can be easily controlled based on their relative overall flow resistance. If the channel dimensions are the same within the structure, the reservoir filling sequence depends on the length of the channel connecting the reservoir and the entry port.

FIG. 2 illustrates another exemplary embodiment in which the reservoirs 8 are connected in series or branched-series configuration. This configuration is useful for applications in which it is desired to add a reagent at different dosages to the sample stream before it reaches different branched reservoirs in the branched-type configuration. For example, if an acidic reagent soluble to the sample stream immobilized in reservoirs A, B, C, and D, liquids in reservoirs E, F, G, and H will contain a different amount of acid. The configuration shown in FIG. 2 makes it possible to study a reaction between the liquid sample with a reagent in branched reservoirs F, G, and H under different pH conditions. This type of fluidic manipulation is usually very difficult to achieve with conventional methods based on pumps and valves.

FIG. 3 shows a parallel series configuration. With this configuration, one can immobilize different soluble reagents in reservoirs A, B, C, D, E, and F, and the same reagent in reservoirs a to d. Thus, one can create an array of sensor reactions between the liquid sample with a reagent in reservoirs a, b, c, d, and f in the presence of reagents delivered through reservoirs A, B, C, D, and F.

Reservoirs A to D shown in FIG. 3 can be used to cover the sensor elements while reservoirs a to d is used for sample volume control. By changing the volume of reservoirs a to d, one can control the effective sample delivered to reservoirs A to D.

In a similar fashion to the process described in FIGS. 2 and 3, reservoirs A, B, C, D, E could contain a material or membrane that removes a species or chemical from the liquid flowing in the channels, thus modifying the liquid or removing interferences prior to the liquids arrival at subsequent reservoirs.

FIG. 4 shows a configuration that allows first filling the reaction reservoir 42 with a first liquid sample for a short period, and then withdrawing the first liquid sample to a waste reservoir 48. After the first liquid is withdrawn to the waste reservoir, a second liquid sample can be added to the entry port, and driven to fill the reaction reservoir. In order to achieve these fluidic functions, capillary pressure generated by the waste reservoir 48 should be greater than that by the reaction reservoir 42 and the channel 5 connecting between the reaction reservoir 42 and the entry port 12. The delay time can be controlled by varying the length or/and flow resistance of the delay line 44.

Another important but optional feature of the fluidic delivery device 10 disclosed in this invention is that it is removable from a substrate. The substrate provides the bottom wall of the reservoir, and the sampler binding layer provides a liquid seal to the substrate. This is especially advantageous comparing to many fluidic devices disclosed in the prior art when sensor or reaction elements have to be included in the reservoir. Using this liquid delivery device, the substrate can be prepared independently.

Since the device disclosed in this invention can be separated from the substrate after its use, it can be a reusable device, although it is also suitable to being treated as a disposable component. Also, the substrate can be reusable if analyte-responsive sensor are reversible or regeneratable.

Reference materials can be used to normalize the sensor response. This can be any stable material whose spectrometric properties are not affected by the environmental or system parameters experienced by the array, e.g., temperature, light, and humidity. Alternatively, they could be the substrate itself whereupon sensor elements are deposited, or these reference materials could be incorporated into the films, affixed to array structure, or be the array material of construction. These materials can be any spectral standard from black to white, of any applicable wavelength suitable to the particular array system design. The reference materials could also be dyes, organic, or inorganic pigments that have spectral bands that do not significantly overlap with spectral bands of the sensor element. The reference material may also comprise an optical response material such as inorganic, organic, and polymeric photonic crystals.

Normalization using the response from the internal references is useful for reduction of errors caused by variation in the optical path length, dimensions of the sensor element, and other variation sources known in the art. More specifically, including an internal reference in the sensor element is important in two situations. First, if the sensor element before exposure is transparently colorless, optical measurements before exposure to the sample may not provide any useful information for sensor response normalization. Secondly, if a measurement on the sensor element cannot be carried out before its exposure to the sample, reading at the λmax (maximum wavelength) of the internal reference after the exposure can be used to correct the sensor response at the λmax of the sensor element. Since the internal reference is an integrated part of the sensor element, tracking changes in its optical response provides information on the physical state of the sensor element after its exposure to the sample and environment. For example, changes in the physical state of the sensor element due to swelling or loss of transparency contribute to the overall sensor response. Differences in signal readings at the λmax of the internal reference before and after the exposure can be used to separate the sensor response due to analyte-sensor interaction from that due to changes in the sensor physical state.

Multiple reference materials can also be deposited on the sensor array and RGB values measured from the reference areas can be used to normalize the sensor response and eliminate any variation that may be caused by illumination changes during the image capture process. The normalization can reduce array-to-array variations introduced in the manufacturing, storage, or sample application processes.

Secondary effects can limit the performance of an array detection system. These effects include noise from the array sensor system, effects from environmental or system parameters, defects caused during the manufacturing or sample application processes, as well as unexplained outliers in the data set, such as interferences, that alter the true analyte response. Minimizing secondary effects can be accomplished by using an individual reduction tool or by combining more than one tool, if appropriate.

Noise reduction can be used to improve array response using several categories of data manipulation. In one envisioned form of noise reduction a digital image file is generated and stored in a computer or microprocessor, and noise reduction methods are applied to analyze the raw data. These can include, but are not limited to, Fourier transforms, wavelet transform, Kalman filtering, Savitsky-Golay smoothing, running mean, median, and polynomial methods. In the case of color response, RGB values over each sensor region can also be averaged. In another instance, selective data elimination can be applied where one calculates standard deviations within a smaller area centered at each pixel inside the sensor element, referred to as subset standard deviation. If the subset standard deviation of a group of pixels is greater than a pre-set value, this group of pixels can be rejected from the set.

In another envisioned form of secondary effect reduction, the system can have elements that sense environmental variables such as temperature, where the temperature measure can be used to account for predetermined variations due to this or similar environmental measures. Similarly, additional measurement elements can be included to account for general system parameters such as sample clarity, system conductivity, oxidation reduction potential, or similar variables that can affect array response. Having measures of these additional variables allows the system to compensate for effects not fully compensated using simpler analytical tools.

In another envisioned form of secondary effect reduction, a system can be established to eliminate response from defective solid sensor film defects. Some defects are caused by out of spec formulations that are used to prepare the sensor films that result in spatial inhomogeneities of the films, or might be introduced into the film preparation steps, such as inclusion of dust particles in the film. Foreign materials might also be deposited on the film during the exposure to the sample matrix. A digital image provides very high spatial resolution on color intensity distribution over each sensor region. This spatial information can be exploited for noise reduction. A variety of data analysis tools can be used to reduce errors from film defects, and an algorithm may be applied to discriminate noises caused by the defects. For example, averages and standard deviations of the RGB values for the whole area of the sensor element are first calculated. We refer to them as set averages and standard deviations, respectively. Then the RGB values from each pixel are compared with the set averages. If the difference is greater than a pre-set multiple of the set deviation, this pixel can be rejected from the set. A similar calculation can be used to reject a group of pixels. Defects in the sensor elements may also exhibit unique color or/and spatial patterns, such as lines and dots. Pattern recognition algorithm may be applied to identify the defect regions. Additionally, defects in the sensor elements are not normally distributed. The optical response from the defect regions is either greater or less than the set average. Thus, normality test can also be used to reject readings from the defects. This is especially useful when the overall sensor film quality is poor and a significant amount of white noise exists in the set.

Interacting system contaminants can also cause errors in array results. Compensation for interferences can be made if the concentration of the interfering species can be measured directly, or if it can be inferred from parallel response from dissimilar sensor films. The interdependent nature of chemical species in solution could result from interferences, where these interferences can be caused by competing reactions of interfering species with the sensing reagent. The traditional wisdom has been focused on the development of interference free chemical reagents for an individual analyte. Chemometric data analysis algorithms have been used to analyze overlapping spectral response for interference reduction, and this has been described in the literature.

A three-part data generation and analysis method is used in this invention to solve interference problems. First, we design sensors to measure parameters that define the chemical and physical state of the sample. These parameters include temperature, pH, and alkalinity. Second, we design sensors that respond independently to a group of interference species. Third, we design sensor films that respond to the same analyte, but that have dissimilar interference response. The sensor responses from these sensors are deconvoluted to review true concentration for each analyte among the interference species. One can also compare the response pattern of the sensor suite from the measured sample to the stored model. The stored model is built from the responses of the sensor films to a range of analyte species and their combinations with the additional responses of the sensor films to expected interferences at their different levels. By capturing different sensor response at various combinations of analytes and interferences, the model captures the response pattern over the analytes dynamic range of interest. The tools for quantitative analysis of sensor films in their combination include neural network, principal components regression, locally weighted regression, partial least squares and any others known in the art.

Multivariate analysis has been widely used in analytical chemistry, especially in spectroscopic analysis. One aspect of the present invention is that a systematic method is utilized to simultaneously determine multiple analyte concentrations in a water or process sample. An exemplary method disclosed in this invention provides a sensor array comprising multiple sensor elements that are chosen to de-convolute the inter-dependent nature of the chemical equilibrium in water or process systems. The sensor array provided may include sensor elements that are specifically designed to be responsive to water or process parameters that are needed for multivariate analysis, and that would otherwise not be needed as part of a simple, but less accurate analysis. Additionally, the array platform allows systems with multi-response chemistries to be used, and that are interpreted by deconvoluting the dual response results.

Determination of pH, alkalinity, hardness, and phosphate are non-limiting examples of the complex nature of a water system testing, where different analytes produce interdependent responses. pH is by definition a measure of hydrogen ion activity that is defined by thermodynamic properties of the water sample. pH is also affected by carbonate concentrations in the same water sample, and carbonate exists in several aqueous forms whose proportions are determined by a complex series of equilibria as defined by the system pH. Carbonate and corresponding phosphate equilibria provide buffering environments. Buffer is a mixture of a conjugate acid-base pair that can resist changes in pH when small amounts of strong acids or bases are added. The buffer capacity of a solution is the number of moles of strong acid or strong base needed to change the pH of 1 Liter of buffer solution by 1 pH unit. Hardness is referred to as the total calcium and magnesium concentration, including numerous forms of calcium and magnesium species that may exist in the system. Some of these calcium and magnesium forms may include phosphates, and these phosphate salts are in equilibrium with soluble forms of the contributing ions. The ion concentrations exist in a complex equilibrium that balances the pH with carbonate, phosphate, and hardness concentrations. Phosphate can also exist in additional forms in water, and again, the respective phosphate forms are determined by a pH, alkalinity, and counterion balance in a series of interrelated equilibria. One can use a phosphate optode to measure phosphate, but a phosphate optode may respond to just one ionic form of phosphate in water, e.g., mono-hydrogen phosphate ions (HPO₄ ²⁻) species. In order to obtain the total phosphate concentration, the sample pH, carbonate, and hardness concentrations must also be known. All the species in the water system are in chemical equilibrium, and all the contributing equilibria must be accounted for when determining a single analyte concentration. One must measure all these analyte concentrations as well as account for environmental properties like temperature to make an accurate measure of the single analyte.

The mathematical details of this complex chemical and thermodynamic equlibria are well known in the art and are beyond the scope of this disclosure, but a simpler example based on alkalinity and pH will give a clearer account of the difficulties associated with multi-equlibria measurements made with solid film sensor that use color responsive reagents. The purpose of this non-limiting example is to illustrate one of the systematic methods disclosed in this invention.

pH is defined by equation below: pH=−log 10a _(H+), (where a_(H+) is activity of hydrogen ion). Hydrogen ion is related to other chemical species in the system through the following equilibrium: H₂O(l)

H⁺(aq)+OH⁻(aq), HA(aq)

H⁺(aq)+A⁻(aq),

HA(aq) stands for a aqueous Bronsted acid, and A⁻ is the conjugate base of HA(aq). The existence of Bronsted acids and bases give rise not only to acidity or basicity of the system but also to pH buffer capacity. The pH buffer capacity is usually measured as alkalinity in the water treatment industry, and this is primarily a function of the total carbonate concentration.

When a sample is applied to a pH sensor region, a pH sensing reagent such as a pH indicator dye, referred to as ‘Ind’ hereafter, interacts with the hydrogen ion in the sample through equilibrium below: Ind+H⁺(aq)

IndH⁺.

The change in the indicator concentration according to above equilibrium is used to determine the pH value of the sample. The indicator molecules (Ind and IndH⁺) have different spectrum, and a change in spectral absorbance indicates a shift in equilibrium, which can reflect shifts in system pH, and the change of the indicator concentration is usually measured by a change in the optical properties of the sensor region. The optical properties include absorbance and fluorescence.

Since the pH indicator itself is usually a Bronsted acid or base, as indicated in the above equilibrium, the pH measuring process perturbs the acid-base equilibrium in the sample, which results in a measurement error for pH. The numerical value of this error is a function of the buffer capacity of the system. Therefore, one has to know alkalinity of the system to determine pH accurately.

A sensor array for pH and alkalinity analysis comprises multiple sensor elements. Some elements measure sample alkalinity, while other elements measure sample pH, and the combination of these multiple sensor elements are used to extend the detection range for the array. The response of the sensor element to alkalinity can be made independent of the sample pH, and the alkalinity of the sample can be obtained from just alkalinity elements. As mention above, the response of the pH sensor elements is a function of both pH and alkalinity. A two-dimensional calibration surface can be obtained for pH and alkalinity. The pH value of the sample can be determined using the measured alkalinity value and interpreting the sensor response using the two-dimensional calibration surface.

Kinetic Response

Often, for quantification, a sensor response should reach a steady state upon exposure to the sample. In practice, some sensors have a long response time and it takes an unacceptably long time to reach steady state, and measuring the sensor response at any single time may result in errors due to variations in timing. For a sensor array, different response reading methods should be applied. For non-steady state sensors, time-dependent measurements are required. Kinetic information can be interpreted for the dynamic characteristics of the system such as the initial slope, slope at a given time, and intercepts of a selected segment of the response curve.

Additionally, temporal response can also provide a sensitive measure of both analyte response as well as reflect the presence or concentration of contaminants that effect the sensor response kinetics. Similarly, the kinetic response can be used to measure the concentration of catalyzing agents that may be in the system, making the temporal response independent of the sensor array equilibrium measure. Many time series statistic models can be used to treat the non-steady state sensor response. In general, the final reading and the responses prior to the final reading can be fitted into a model to minimize the instrumental and measurement errors.

Total Analysis System

The total analysis system of the present invention includes an optical array platform comprising diverse chemically or physically responsive sensor films. The system produces an optical response proportional to the desired chemical or physical parameter, provides secondary effect reduction from noise, defect, and interference effects, compensates for multivariate interactions, accounts for test array history, and provides a reference system to calibrate the sensor array response to the optical detection platform. This complex test array can be combined with time-based data acquisition to provide temporal test analysis that can further enhance overall array response. The complex array elements described herein show how each element enhances the array performance, and how the combination of these elements produces optimized environmental and biological measurements. Further, this enhanced optical array platform is advantageously suited for non-laboratory environments.

In other aspects of the invention, internal reference materials are included as an integrated part of the sensor array. The reference materials allow normalization of sensor response to eliminate variations caused by variation in illumination, sensor element quality, and environmental parameters. In addition, the present invention provides solutions to specific problems associated with the digital imaging techniques known in the art.

In other aspects of the invention, the sample volume that each sensor element is exposed to is controlled by a capillary-flow-based fluidic device which transports and doses, in metered quantities, a controlled volume of sample liquid to the sensor elements. In this way, the fluid delivery device of the present invention makes it possible to efficiently construct a sensor array with both reversible and irreversible sensor elements.

The following examples are included to demonstrate the broad applicability of the present invention. It should be appreciated by those of skill in the art that the techniques disclosed in the examples which follow represent techniques discovered by the inventors, and thus can be considered to constitute exemplary modes for its practice. However, those of skill in the art should, in light of the present disclosure, appreciate that many changes can be made in the specific embodiments which are disclosed and still obtain a like or similar result without departing from the scope of the invention.

EXAMPLE 1

Reservoir fill time as a function of channel and reservoir geometric parameters. The channel and reservoir layout are shown in FIG. 20. The device comprises three layers as shown similarly in FIG. 1. The top cover layer 2 is a heat sealable hydrophilic film. Vent holes 7 (1.5 mm diameter) were cut through this layer. The middle channel layer 4 is a 0.78 mm thick polycarbonate sheet with open channels 5 and rectangle openings (i.e., reservoirs) 8 cut by means of computer numerical controlled (CNC) machining. The bottom sampler binding layer 6 is a 40 Shore A silicone gasket, providing seal to a substrate. Rectangle openings are die cut through the gasket. When these layers are laminated to form a fluid delivery device, channels are created between the top hydrophilic layer 2 and middle channel layer 4. The rectangle openings of the channel layer and sampler binding layer define an open-bottom reservoir, with the hydrophilic layer as its top wall. When this assembly is attached to a substrate, closed reservoirs are formed and connected to a central sample entry port through channel.

Channel and reservoir parameters tested in this example are listed below in Table 1. It was found that reservoir fill time (t/seconds) can be expressed a function below: Log(t)=K−0.97251 Log(W)−2.43118 Log(D)+1.34630 log(L)+1.70630*Log(D_(gaset))  Equation (1) where L, W, and D are channel length, width, and depth, respectively; D_(gasket) is the gasket thickness. Constant K equals −1.9944 for single channel to single reservoir configuration, and −1.7740 for single channel to two reservoirs configuration. FIG. 21 shows the comparison of fill time predicted by the equation above with experimental data.

TABLE 1 Ranges of channel and reservoir parameters studies in Example 1. Reservoir per W/mm D/mm L/mm Gasket/mm channel 1 0.2 45 0.47 1, 2 1 0.3 45 0.47 1, 2 2 0.2 45 0.47 1, 2 2 0.3 45 0.47 1, 2 1 0.2 45 0.55 1, 2 1 0.3 45 0.55 1, 2 2 0.2 45 0.55 1, 2 2 0.3 45 0.55 1, 2 1 0.2 12 0.55 1, 2 1 0.3 12 0.55 1, 2 2 0.25 12 0.55 1 2 0.25 21 0.55 1 2 0.25 30 0.55 1 2 0.25 38.5 0.55 1 Reservoir depth = Gasket + 0.78 mm Reservoir length = 6.0 mm Reservoir width = 5.0 mm

Based on this equation, one can design a device to allow all the reservoirs being filled in a narrow time range although the distance of a reservoir to the central entry port varies. If it is desirable, reservoirs can be filled sequentially by choosing channel parameters according to the equation above.

EXAMPLE 2

A 54-reservoir sample delivery device is shown in FIG. 1. The device 10 is assembled from four components by a method similar to that described in Example 1. The channel 5 depth and width are 0.33 mm and 1.5 mm, respectively. The length and width of reservoirs 8 are 5.5 and 4 mm, respectively. The sampler binding layer 6 is cut from a 0.55 mm thick, adhesive backed, and transparent silicone rubber sheet. The thickness of the polycarbonate channel layer 4 is 0.78 mm. Selection of these design parameters was guided by the equation shown in Example 1. A 2.7 ml sample solution containing 100 ppm Basic Blue was delivered to the sample entry port. Real time flow in the channels and reservoirs were monitored using a digital video camera. Fill time for each reservoir was retrieved from the recorded video films. Average fill times for all 54 reservoirs obtained from six devices are presented in FIG. 22. The data demonstrate that the device allows delivering a liquid sample to multiple reservoirs in a narrow time range.

EXAMPLE 3 Determination of Chlorine Concentration in a Water Sample

Six chlorine sensitive films were deposited on a thin translucent polyethylene sheet. A 20 μl chlorine standard solution, prepared from a 5% NaOCl by dilution with deionized water, was spotted on each film. The water sample was removed from the films 1 minute after spotting. A blue color was developed as chlorine in the water sample reacts with the chlorine sensitive reagent immobilized in the films. The image of these six films was captured with Hewlett Packard scanner ScanJet 6300C and is shown in FIG. 5. The digital file produced from the scanner was in the JPEG format (67KB). The color depth was 255. The pixel resolution is 200 dpi.

The digital image was processed with Adobe Photoshop® 6. The film areas were selected using the selection tools provided by the Photoshop software package. Average RGB values for each selected color area are listed in Table 2 below. RGB values for the white paper area of the image, referred to as R_(w), G_(w), and B_(w), are also listed in Table 2.

As shown in FIG. 6, the quantity defined in Equation 2 below can be used to quantify chlorine concentration. R_(chlorine)=−log(R/R_(w))−log(G/G_(w))−log(B/B_(w))  Equation (2)

TABLE 2 Color analysis of chlorine sensitive films shown in FIG. 5 Log(R_(w)/R) + Chlorine/ log(G_(w)/G) + ppm R G B log(B_(w)/B) 1 233.63 242.67 246.35 0.0724 2 228.52 243.90 241.61 0.0883 4 219.94 242.00 236.37 0.1178 5 211.46 238.44 235.95 0.1421 10 184.87 228.44 230.34 0.2295 50 124.76 199.12 190.77 0.5418 White paper 254.56 254.60 254.62 0.000 background

EXAMPLE 4 Determination of Alkalinity with Multiple Sensor Regions

Six alkalinity sensitive films were deposited on a glass slide. Suitable sensor types used for this example are described in our co-pending patent applications entitled “Material Compositions for Sensors for Determination of Chemical Species at Trace Concentrations and Method of Using Sensors” and “Self-Contained Phosphate Sensors and Method for Using Same” filed on the same date as the present application, and will not be repeated herein. Unlike chlorine analysis, multiple films were used for determination of alkalinity of a single water sample. A 20 μl alkalinity standard solution was spotted on each of the six films. The water sample was removed from the films 2 minutes after spotting. Ten alkalinity solutions were measured.

As shown in FIG. 7, the image of total 60 exposed films was captured with Hewlett Packard scanner ScanJet 6300C. The digital file produced from the scanner was in the JPEG format (48KB). The color depth was 255. Absorbance of the each exposed film at 650 nm was measured with an Ocean Optics USB2000 spectrophotometer.

Average RGB values for the selected color areas are listed below in Table 3. RGB values of the white paper background are 239.41, 239.34, and 244.19, respectively. The following quantity Ralk is used to quantify alkalinity: R_(alk)=[(R_(w)−R)²+(G_(w)−G)²+B²]^(1/2)/B_(w).  Equation (3)

TABLE 3 Color analysis of films shown in FIG. 7 Alkalinity/ppm Film number R G B R(alk) 18.8^(a) 1 227.44 227.46 137.72 0.5682 2 226.92 226.73 140.2 0.578725 3 221.63 225.69 135.09 0.560781 4 221.69 225.33 135.85 0.563968 5 219.89 224.98 136.94 0.569506 6 216.71 225.18 143.37 0.59726 Average 0.573073 37.7^(a) 1 221.76 225.26 139.83 0.580045 2 194.52 218.29 168.23 0.718228 3 202.97 220.22 152.49 0.646812 4 195.72 219.17 165.31 0.705072 5 187.68 215.39 157.51 0.685975 6 179.48 214.16 162.75 0.717686 Average 0.675636 75.5^(a) 1 180.09 212.77 157.7 0.698513 2 170.68 211.28 175.25 0.779415 3 159.29 206.99 178.37 0.811646 4 152.1 204.67 185.54 0.851659 5 149.23 203.04 186.75 0.862184 6 152.45 203.97 187.63 0.859187 Average 0.810434 39.0^(b) 1 210.27 220.64 132.93 0.562534 2 192.88 216.61 150.25 0.65082 3 189.62 216.16 159.76 0.691825 4 182.28 213.28 153.93 0.680802 5 176.07 212.64 169.23 0.748012 6 183.77 214.95 173.86 0.754201 Average 0.681366 58.5^(b) 1 187.5 215.83 154.64 0.674907 2 176.74 213.27 178.27 0.781173 3 174.38 212.25 172.24 0.762069 4 167.33 210.12 176.22 0.788816 5 164.12 208.65 178.48 0.803171 6 171 210.72 193.42 0.848307 Average 0.776407 97.5^(b) 1 158.03 206.39 193.39 0.86976 2 150.32 203.24 188.28 0.865716 3 138.93 198.84 189.03 0.892229 4 139.27 198.77 192.7 0.90472 5 135.26 196.68 194.44 0.920039 6 133 195.89 198.3 0.938624 Average 0.898515 117.0^(b)  1 142.43 200.58 197.08 0.913398 2 139.2 198.84 193.74 0.908514 3 133.99 196.92 198.26 0.935815 4 132.53 195.84 195.89 0.931041 5 131.33 195.88 199.17 0.9449 6 132.71 196.04 196.04 0.931068 Average 0.927456 156.0^(b)  1 132.22 195.53 197.5 0.937564 2 127.95 193.56 202.14 0.963713 3 123.56 191.45 200.69 0.969018 4 110.08 189.51 205.3 1.014391 5 119.48 189.62 202.06 0.983554 6 124.03 192.11 203.77 0.97827 Average 0.974418 60.0^(c) 1 187.97 215.12 143.39 0.631684 2 168.96 209.56 165.01 0.744808 3 160 207.17 168.27 0.773279 4 159.19 207.64 180.98 0.821017 5 155.67 206.24 178.99 0.820522 6 158.72 208.13 192.23 0.863269 Average 0.775763 121.9^(c)  1 148.73 203.1 195.82 0.896101 2 140.39 199.41 195.24 0.911284 3 138.33 198.51 198.69 0.928096 4 132.01 195.96 199 0.942936 5 130.44 194.99 199.12 0.947129 6 133.34 196.93 201.73 0.949377 Average 0.929154 ^(a)Solutions prepared from Na₂CO₃ ^(b)Solutions prepared from a mixture of Na₂CO₃, NaHCO₃, and Na₂HPO₄. The percentage alkalinity contributions of Na₂CO₃, NaHCO₃, and Na₂HPO₄. are 10%, 80%, and 10%, respectively. ^(c)Solutions prepared from NaHCO₃

In FIG. 8, the average R_(alk) value of the six films is plotted as a function of solution alkalinity. Note that the calibration curve for alkalinity analysis is not necessarily a straight line. The curvature in the calibration curve is not due to the current color analysis method. This is supported by the linear correlation of R_(alk) to absorbance measured with the spectrophotometer at 650 nm as shown in FIG. 9.

EXAMPLE 5 Determination of Alkalinity Using a Digital Camera—Normalization Against Internal Reference Areas

The image of total 60 films was captured with Sony DSC S75 digital camera. The camera was set in the automatic mode in which white balance, focusing, and aperture were automatically adjusted. The glass slides were placed near a 40-Watt desk lamp.

Average RGB values for the selected color areas are listed below in Table 4. Unlike the image captured by a digital scanner, illumination cross the subject is not uniform. In order to correct this, RGB values of a white paper background near each film were taken. Instead using a single set of RGB values for the white background in Equation 3, each color film has a set of R_(w)G_(w)B_(w) values, as listed in Table 4 below. The calibration curve is shown FIG. 10.

TABLE 4 Color analysis of the digital camera image for Example 5 Film # R G B Rw Gw Bw R(alk) 18.8 ppm 1 140.86 143.1 84.25 162.69 157.08 135.4 0.651018 2 140.67 142.3 82.65 165.08 159.68 137 0.641711 3 146.04 147.05 84.465 168.13 161.59 140.64 0.629325 4 147.2 147.73 85.82 170.18 164.35 142.98 0.632149 5 153.32 153.27 89.03 172.23 166.67 144.92 0.634814 6 154.87 156.16 88.95 174.44 169.12 147.71 0.622807 Average 0.635304 37.7 ppm 1 131.98 146.95 98.97 169.39 164.92 144.31 0.743674 2 134.93 144.73 95.6 168.45 162.99 142.68 0.721466 3 136.82 144 89.99 167.09 161.07 141.07 0.683823 4 137.45 141.81 88.41 164.71 158.89 137.82 0.682634 5 133 140.36 96.32 162.02 156.08 134.82 0.755211 6 137.23 138.96 81.32 158.58 153.87 131.44 0.649633 Average 0.706073 75.5 ppm 1 116.72 144.57 110.98 174.52 168.4 146.2 0.871262 2 117.94 148.16 114.66 177.59 170.89 149 0.880748 3 120.48 150.52 115.13 178.24 172.04 150.85 0.865707 4 123.11 153.27 113.07 180.76 174.86 153.65 0.837891 5 132.76 156.82 111.64 182.99 177.51 156.33 0.794189 6 140 159.85 107.86 184.27 179.33 158.15 0.747441 Average 0.832873 39.0 ppm 1 157.33 167.69 116.34 187.75 181.8 161.13 0.75142 2 152.09 166.7 113.08 186.53 180.76 159.11 0.748171 3 150.41 163.15 109.64 185.09 179.19 157.94 0.735136 4 152.52 161.57 107.71 182.21 176.47 154.05 0.731686 5 153.53 160.2 107.31 180.93 174.54 152.83 0.730729 6 154.42 158.39 96.93 178.45 172.23 150.17 0.671364 Average 0.728084 58.5 ppm 1 130.39 144.37 108.43 169.76 163.45 141.56 0.825964 2 130.02 147.89 110.69 173.18 166.58 145.38 0.827266 3 125.14 148.35 105.91 174.06 168.29 146.3 0.808982 4 133.16 152.51 105.58 176.9 170.27 148.63 0.778131 5 137.64 155.4 111.67 178.51 172.87 152.09 0.790259 6 148.67 158.85 102.81 179.85 175.19 154.16 0.704914 Average 0.789253 97.5 ppm 1 106.9 154.16 124.43 181.09 175.61 154.76 0.946293 2 107.42 152.28 120.38 179.55 174.12 152.9 0.928874 3 109.2 150.24 117.63 177.87 171.84 150.24 0.917926 4 107.62 146.64 113.12 174.96 169.24 147.67 0.904533 5 124.06 148.49 113.67 174 167.34 147.23 0.852947 6 118.91 143.63 110.02 170.08 163.88 142.39 0.863934 Average 0.902418 117.0 ppm 1 102.02 139.66 111.68 170.9 164.47 143.33 0.931683 2 99.01 141.84 112.46 174.03 167.55 146.01 0.942463 3 103.35 144.47 114 175.36 168.9 147.89 0.926593 4 106.92 147.08 117.41 177.91 172.32 149.94 0.930408 5 115.69 152.03 117.52 179.77 174.24 152.46 0.889973 6 121.83 155.64 120.84 181.74 176.32 155.59 0.876998 Average 0.916353 156.0 ppm 1 114.04 147.2 117.91 174.8 168.63 146.92 0.914541 2 104.21 148.52 121.67 178.71 173.28 150.76 0.960464 3 100.84 150.88 121.22 180.18 174.49 152.95 0.959709 4 100.89 153.52 126 182.97 177.44 155.97 0.976259 5 101.37 155.15 127.03 184.47 179.49 158.37 0.970738 6 104.23 157.53 130.45 186.1 180.74 160.09 0.972901 Average 0.959102 60.0 ppm 1 122.73 135.27 98.11 160.28 154.92 131.5 0.812717 2 116.3 136.7 94.48 163.61 157.33 135.23 0.796112 3 120.64 139.51 99.26 165.64 158.85 138.02 0.801963 4 127.22 142.71 95.65 168.63 161.2 141.11 0.75017 5 133.11 146.56 97.41 169.92 163.69 141.93 0.743554 6 144.59 151.32 93.69 170.84 165.67 144.44 0.680908 Average 0.764237 121.9 ppm 1 115.97 158.69 127.09 176.25 170.76 149.63 0.943514 2 113.09 157.26 124.81 172.85 169.22 148.45 0.935635 3 115.71 155.67 121.02 174.25 168.47 147.28 0.916913 4 117.18 154.2 120.59 172.42 166.52 144.39 0.922578 5 117.79 151.55 117.27 170.7 163.44 141.64 0.912184 6 122.59 149.93 113.83 167.05 159.76 138.39 0.885897 Average 0.919454

The image in this example was captured about 36 hours after the image in Example 2 was taken. In order to compare the digital camera results with those obtained with scanners, another image was prepared using Canon N650U scanner at the same time when the camera image was taken. Correlation of the relative blue intensity from the Canon scanner to that from Sony digital camera is shown in FIG. 11. The result shown in FIG. 11 indicates that normalization against an internal color reference can effectively eliminate errors caused by non-uniform illumination when a camera is used to capture the sensor image.

EXAMPLE 6 pH Determination by Means of Multivariate Calibration

A pH sensitive film was deposited on a polycarbonate sheet. The film contains pH indicator dye bromothymol Blue and other additives. The pH standard solutions used in this example were prepared from sodium carbonate and sulfuric acid solutions. A glass electrode, calibrated against two pH buffers (7.00 and 10.00 from Fisher Scientific, traceable to NIST standards), was used to measure pH values for the standard solutions. No correction was made for the ionic strength effects on hydrogen ion activity coefficient and liquid junction potentials. Alkalinity was measured by titration against a 0.2 N sulfuric acid solution.

A 40 μl aliquot of the pH standard solution was spotted on the film. The sample was removed after 2 minutes, and the spotted area was dried by a modest blow of airflow. Digital images of the film before exposure and after exposure were captured with using Canon LiDE 80 scanner in 48-bit color mode with 300 dpi spatial resolution. The image file was saved in the uncompressed TIFF format. Photoshop CS was used to retrieve RGB values from the file. The RGB values listed below in Table 5 were averaged over a 1000-pixel square centered at the spotted area. The following quantity R_(pH) is chosen as a sensor response to quantify sample pH: R_(pH)=(R/G−B/G)_(exposed)−(R/G−B/G)_(unexposed)  Equation (4)

TABLE 5 Color analysis of digital images for Example 6 Unexposed Exposed Response PH Alkalinity/ppm R G B R G B R_(pH) 8.34 97.4 93.45 175.28 158.23 87.61 157.64 139.94 0.3705 8.48 100.6 94.67 177.02 159.89 109.57 161.73 129.58 0.3695 9.07 112.4 94.06 176.91 159.30 117.92 164.81 128.52 0.3365 9.36 92.4 92.97 175.08 158.03 120.91 166.85 127.96 0.2917 9.43 103.2 93.95 176.48 159.21 136.19 173.50 128.34 0.2778 10.18 100.0 93.67 176.33 158.83 125.32 165.97 125.32 0.1106 9.80 303.2 93.84 176.33 158.82 124.92 169.49 130.35 0.0376 8.94 311.0 93.23 175.63 158.26 116.49 165.61 129.50 0.2447 8.40 303.6 94.13 177.03 159.34 115.52 165.44 130.51 0.3045 7.97 301.8 102.17 178.73 158.17 101.40 161.05 134.05 0.3293 7.08 308.6 93.67 176.33 158.83 125.32 165.97 125.32 0.4150 7.04 510.8 94.91 176.01 156.53 118.66 156.56 116.85 0.3617 8.18 516.8 95.42 176.31 156.27 100.98 155.56 125.36 0.1884 8.54 506.0 98.57 177.33 156.99 102.86 153.85 123.80 0.1933 8.84 490.0 94.21 175.45 155.48 97.51 154.32 126.18 0.1634 9.12 504.0 99.25 177.38 157.83 95.49 158.61 135.89 0.0755 9.48 500.0 101.78 177.08 155.51 86.71 150.73 133.38 −0.0062 9.63 492.0 97.58 176.17 157.58 82.32 151.37 138.21 −0.0286

The sensor responses are plotted as a function of sample pH and alkalinity in FIG. 12. It is clear from FIG. 12 that the sensor response is a function of both sample pH and alkalinity as described in the above section. It was found that experimental pH values could fit into the following two-variable calibration equation within 0.09 pH unit (average absolute deviation). pH=a₀+a₁alk+(a₂+a₃alk) R_(pH)+(a₄+a₅alk)(R_(pH))²  Equation(5) Values of fitting parameters a0 to a5 are listed below in Table 6. The pH values calculated from the above equation are compared with experimental values in Table 6.

TABLE 6 Parameters for the pH calibration equation and pH values calculated from the calibration equation Alkalinity/ Response ppm pH_calculated pH_Measured Difference 0.3705 97.4 8.45 8.34 0.11 0.3695 100.6 8.45 8.48 0.03 0.3365 112.4 8.83 9.07 0.24 0.2917 92.4 9.39 9.36 0.03 0.2778 103.2 9.47 9.43 0.04 0.1106 100.0 10.25 10.18 0.07 0.0376 303.2 9.73 9.80 0.07 0.2447 311.0 8.86 8.94 0.08 0.3045 303.6 8.38 8.40 0.02 0.3293 301.8 8.14 7.97 0.17 0.4150 308.6 7.10 7.08 0.02 0.3617 510.8 6.98 7.04 0.06 0.1884 516.8 8.42 8.18 0.24 0.1933 506.0 8.43 8.54 0.11 0.1634 490.0 8.70 8.84 0.14 0.0755 504.0 9.16 9.12 0.04 −0.0062 500.0 9.52 9.48 0.04 −0.0286 492.0 9.60 9.63 0.03 Parameters for equation 3 in example 4 A0 A1 A2 A3 A4 A5 10.067 −0.00113 7.8279 −0.02317 −20.2109 0.04331

EXAMPLE 7 Noise Reduction

Several solid sensor films have been produced with defects that were caused by undissolved reagents of the sensor film formulation prepared from a polymer solution. Some defects were introduced in the film preparation steps, such as inclusion of dust particles in the film. Foreign materials could be deposited on the film during the exposure to the sample matrix. A digital image provides very high spatial resolution on color intensity distribution over each sensor region. This spatial information can be exploited for noise reduction. A variety of data analysis tools can be used to reduce errors resulted in by the film defects. This example demonstrates a simple statistical approach to reject readings from the defect areas.

An enlarged color image of a pH sensor spot taken from Example 6 is shown in FIG. 13. A dust particle, indicated by the reference number 130, is visible in this image. RGB values for a 40-pixel horizontal line containing this dust particle are shown in FIG. 14. In Example 6, RGB values averaged over the whole sensor region are used to calculate the sensor response. For simplicity, we use one-dimensional data in this example to demonstrate the noise reduction method. We first calculate averages and standard deviations of this data set respectively for R, G, and B. Then, we reject those points with deviation from the set average greater than a pre-set multiple of the standard deviation for the set. Finally, averages and standard deviation are calculated from pixels excluding the dust areas.

A similar calculation can be used to reject a group of pixels. We first calculate the standard deviations for the whole sensor regions. Then, we calculate standard deviations within a smaller area (a 6-pixel circle chosen for this example) centered at each pixel along each horizontal pixel line. FIG. 15 shows results from these calculations. It is clear that an area centered at the 40th pixel should be rejected.

FIG. 16 shows a calibration curve for a calcium sensor film. The film was prepared from a polymer solution, which contains a calcium responsive dye. The film was prepared on a polycarbonate sheet with a film applicator. When the film is dry, some dye aggregate to form small dark areas randomly distributed across the entire film, barely seen by naked eyes. A digital image of the exposed film was scanned with Canon LiPE 80 scanner in 16-bit color mode with 300 pdi spatial resolution. The same equation as used for pH calibration was used to calculate the sensor response. The data filtering methods described above (2× standard deviation) was used to rejected data points originated from the dye aggregates. Generally, 90 to 145 pixels are rejected among a 2700 pixel area. The R-squared value of the calibration is 0.9930, which is significantly improved compared to 0.9886 obtained with un-filtered RGB values retrieved using Photoshop CS.

EXAMPLE 8 Kinetic Response

A polymer solution containing molybdate sensing reagents was prepared. The polymer solution was deposited on a polycarbonate sheet using a film applicator. The polycarbonate substrate was cut into 20×80 mm strips. The bulk of sensor film on each strip was cut and removed, which left only a 6×6 mm sensor spot on the strip. A 200 μm deep, 6.5 mm wide channel was built using a glass slide and double-sided tapes to cover the 6×6 spot to form a fluidic assembly. This assembly, a Canon LiPE 80 scanner, and a 10-ppm molybdate standard solution were placed in a temperature-controlled room. After an equilibration for one hour, sample solution was introduced to the sensor film through the channel by capillary action. Images of the sensor film were acquired at time intervals shown in Table 7 below. Images obtained at 25.3° C. are shown in FIG. 18.

In this example, we want to demonstrate the importance of data analysis for sensor responses, which is an important part of the systematic method for simultaneously determination of multiple analytes disclosed in this invention.

It was found that the initial sensor response of an unexposed film is a function of temperature. For convenience, we normalize the sensor response after exposure by calculating the ratio below to quantify molybdate concentration: R_(Mo)=(R/G−B/G)_(exposed)/(R/G−B/G)_(unexposed)  Equation (6)

TABLE 7 Kinetic data at three temperatures Time/s R G B R_(Mo) T = 4.5° C. 0 198.71 81.40 65.34 1.0000 21 185.96 87.22 84.94 1.4146 68 179.06 88.60 88.85 1.6092 112 173.74 87.29 87.22 1.6530 170 171.18 86.73 86.79 1.6839 221 169.87 87.04 87.22 1.7255 294 168.82 87.03 87.33 1.7498 486 168.50 87.36 88.05 1.7792 612 168.62 87.68 88.51 1.7933 T = 25.3° C. 0 196.73 74.73 56.43 1.0000 40 189.06 71.78 64.43 1.0813 150 154.68 68.59 66.73 1.4642 225 148.57 67.56 65.76 1.5317 302 146.36 67.57 65.83 1.5753 390 144.61 67.31 65.70 1.6014 477 143.61 67.09 65.51 1.6128 575 142.10 67.67 66.06 1.6708 T = 34.5° C. 0 200.69 78.13 56.06 1.0000 34 199.14 76.14 62.32 1.4146 83 189.60 76.14 70.23 1.6092 133 184.29 76.8 74.20 1.6530 221 175.36 74.87 73.91 1.6839 326 165.41 71.69 70.92 1.7255 409 160.84 70.20 69.45 1.7498 514 158.14 69.54 68.75 1.7792 588 157.16 69.49 68.78 1.7933

This quantity is plotted as a function of time in FIG. 17. Unlike many chemical sensor reactions, the sensor response as defined above does not reach a plateau. It rather continues to increase linearly with time. Linear fits of the last four points for each temperature are shown in FIG. 19.

Taking a single reading from this unsteady film response may result in a large error. For this type of sensor response, we intend to use a quantity derived from the kinetic measurement to quantity the analyte concentration. It is found that the intercepts of the linear curves shown in FIG. 17 are a linear function of temperature. Thus, a multivariate-calibration equation, in which temperature and the sensor responses at several exposure times are independent variable, is suitable for this type of sensor.

Those who are familiar with the art would recognize that many statistical and mathematic models could be used to interpret the kinetic data presented in this example. The methods include Kalman filtering, least square fitting, and other time series prediction tools as detailed in analytical literature.

EXAMPLE 9 Sample-volume Controlled Sensor Array

Eight magnesium sensitive sensor films were screen printed on a 127.8×85.0 mm polycarbonate sheet. A sample delivery device, similar to that described in Example 2, was put on the top the polycarbonate sheet to form enclosed channels and reservoirs. The sensor films are 4 mm long, 4 mm wide, and about 0.01 mm thick. The reservoirs are 5.25 mm long, 5.25 mm wide, and 1.6 mm deep. The volume of the reservoirs is 44.1 μl. When sample is introduced to the central entry port, the capillary force drives the sample to fill the reservoirs. This delivery device provides a volume controlled sample distribution means for the sensor array.

A custom-made 4×4 LED/photodiode array detector was used to monitor the sensor film response. The LEDs have emission maximum at 467, 530, and 634 nm. The LEDs and photodiodes were fixed on two separated print circuit boards. The print circuit boards were held parallel inside an enclosure, where the sensor array assembly can be inserted and aligned with the LED/photodiode array.

A 3.0 ml water sample containing 12 to 100 ppm magnesium was first introduced to the sample entry port. The sensor absorbance at 530 nm (G) 634 nm (R) was measured 3 minutes after sample introduction. It was found that the ratio G/R is linear with respective to magnesium concentration in the sample. The calibration curve is shown in FIG. 23.

While the disclosure has been illustrated and described in typical exemplary embodiments, it is not intended to be limited to the details shown, since various modifications and substitutions can be made without departing in any way from the spirit of the present disclosure. As such, further modifications and equivalents of the disclosure herein disclosed may occur to persons skilled in the art using no more than routine experimentation, and all such modifications and equivalents are believed to be within the scope of the disclosure as defined by the following claims. 

1. A method for simultaneously measuring a plurality of chemical or biological analytes in a sample liquid, said method comprising the steps of: providing an optical sensor array comprising a plurality of solid-film sensor elements for responding to at least one of said plurality of analytes; attaching a sample distribution device to the sensor array, the sample distribution device forming a plurality of fluidic channels leading from a sample entry port to a plurality of reservoirs aligned with said sensor elements, wherein a cover layer of the sample distribution device has a hydrophobic outer surface, a hydrophilic lower surface, and a plurality of vent holes aligned with said reservoirs; introducing said sample liquid into said sample distribution device through the sample entry port; driving said sample liquid along said channels by the hydrophilic surface of the cover layer and into said reservoirs; creating a chemical reaction between said at least one analyte and at least one reagent in said sensor elements; filling each reservoir by accumulating the sample liquid over the top hydrophilic surface above the reservoir to form a pendant droplet, allowing gravitational force to enable the droplet to reach a bottom wall of the reservoir, and then using capillary forces originating from the side walls of the reservoir to drive the sample liquid to fill the entire reservoir while air escapes through the vent hole as the reservoir fills with sample liquid, wherein the hydrophobic outer surface prevents sample liquid from flowing through the vent holes such that the vent holes serve as capillary barriers by impeding the sample liquid from flowing through the vent holes to reach the hydrophobic outer surface of the cover layer, thereby allowing a hydrostatic balance to be reached such that a fixed volume of sample liquid is delivered to each reservoir and sensor element without requiring an external force other than gravity to move the sample through the sample distribution device; generating an optical response from said sensor elements; recording said optical response; and utilizing said optical response to measure, monitor, and control concentration of said analytes in said sample.
 2. The method of claim 1, further comprising the step of directing said sample through capillary barriers between said entry port and said reservoirs.
 3. The method of claim 2, further comprising the step of optimizing channel and reservoir geometric parameters to control and meter flow of said sample.
 4. The method of claim 1, wherein said reservoirs are positioned on top of said sensor elements.
 5. The method of claim 1, wherein said channels and reservoirs are configured in a branched-type structure.
 6. The method of claim 1, wherein said channels and reservoirs are configured in a series or parallel-series configuration.
 7. The method of claim 1, further comprising the step of providing a waste reservoir and a delay channel.
 8. The method of claim 1, wherein said solid-film sensor elements are deposited onto a substrate, and said reservoirs are positioned on top of said solid-film sensor elements.
 9. The method of claim 1, wherein said recording step comprises use of a digital or analog imaging device.
 10. The method of claim 1, further comprising the step of processing said optical response by using Fourier transform, wavelet transform, Kalman filtering, Savitsky-Golay smoothing, running mean, median, polynomial methods, or combinations thereof.
 11. The method of claim 1, further comprising reducing interference from said sensor elements by applying multivariate analysis to said plurality of responses from each said sensor element, wherein said reducing step comprises using multiple sensor chemistries for each said analyte, and deconvoluting interference responses from said multiple sensor chemistries using chemometric analysis, response surface analysis, least square analysis or combinations thereof.
 12. The method of claim 1, wherein said array response is dynamic or steady state.
 13. The method of claim 1, wherein said sensor elements are a combination of reversible and irreversible sensor elements.
 14. The method of claim 1, further comprising the step of utilizing a response of one analyte responsive sensor element to enhance analysis of the other analyte-responsive sensor elements.
 15. The method of claim 1, further comprising the step of providing two or more sensor elements responsive to different analytes to determine concentration of common interferants.
 16. The method of claim 1, wherein said optical response is obtained by using reflected light analysis, colorimetric analysis, scattered light analysis transmitted light analysis, fluorescence analysis, chemoluminescence analysis, photoluminescence analysis, polarimetric analysis, or combinations thereof.
 17. The method of claim 16, wherein said optical response is obtained in the spectral range from about 190 nm to about 10,000 nm.
 18. The method of claim 16, wherein said optical response is obtained in the spectral range from about 200 nm to about 2,500 nm.
 19. The method of claim 16, wherein said optical response is obtained in the spectral range from about 200 nm to about 1,500 nm.
 20. The method of claim 1, further comprising the step of repeating steps at predetermined time intervals before, during, and after exposure to said sample to obtain a temporal response, and utilizing said temporal response to determine rates of response for said sensor array, and concentrations of analytes in said sample.
 21. The method of claim 20, further comprising the step of utilizing said temporal response to identify system variables that are reflected in temporal changes in noise or interference information obtained from at least two-sensor array responses to improve control of said concentration of analytes in said substance.
 22. The method of claim 21, wherein said temporal response is used to identify and eliminate outlier responses that do not reflect real system changes.
 23. The method of claim 1, further comprising the step of providing an internal optical reference material to calibrate said array response of said sensor elements.
 24. The method of claim 23, wherein said reference material comprises an optical response material comprising dyes, paints, and/or pigments.
 25. The method of claim 23, wherein said reference material comprises an optical response material comprising inorganic, organic, and/or polymeric photonic crystals.
 26. The method of claim 23, wherein said reference material is incorporated into said sensor elements or is positioned in proximity to said sensor elements.
 27. The method of claim 1 further comprising depositing the plurality of sensor elements on a base plate and attaching the sample distribution device to the base plate, wherein the sample distribution device comprises the cover layer and a fluidic distribution layer, and wherein the fluidic distribution layer has a plurality of grooves and cutouts that form the channels and reservoirs.
 28. The method of claim 27 further comprising providing a double-sided adhesive film with through holes cut to match said reservoirs in said fluidic distribution layer and binding said sensor film carrying base plate with said fluidic distribution layer using said adhesive film to create the plurality of reservoirs, each containing one of said sensor films, and channels connecting said reservoirs to said sample entry port. 